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We report the first observation of B° J/V'/o(980) and first evidence for -5> J/V'/o(1370), 
wliicli are CP eigenstate decay modes. These results are obtained from 121.4 fb~^ of 
data collected at the T(5S) resonance with the Belle detector at the KEKB e^e~ col- 
lider. We measure the branching fractions B{B'^ — >• J/i/'/o(980); /o(980) tt+tt") = 
(l-16lo:ig(stat.)t|];;[^(syst.)t°;^^(iV^(.)^(.))) x 10"" with a significance of 8.4cr, and B{B° 

J/^/o(1370);/o(1370) 7r+^-) = (0.34lO:i^(stat.)«;EJ^2(syst.)tH5(^s(-)s(-) )) x 10"* with a sig- 

nificance of 4.2(t. The last error listed is due to uncertainty in the number of produced _Bs*'_Bi*' 
pairs. 

PACS numbers: 13.25.Hw 



The b — >■ CCS process 5° — ?> J/ip(f>, which has a rel- 
atively large branching fraction, has been used to ex- 
tract the 5° decay width difference AT and CP-violating 
phase from time-dependent analyses [H . The parame- 
ter Ps is expected to be small in the Standard Model and 
can be sensitive to New Physics. The same b ccs pro- 
cess can also produce the decay B'^ — J/ipfo{980), which 
is another promising channel for accessing the mixing pa- 
rameters, with the clear advantage that no angular analy- 
sis is required because of the — 0+ quantum numbers 
of the /o(980). 

Leading-order light-cone QCD predicts that 
the branching fraction B{B", -> J/i/'/o(980)) = 
(3.1 ± 2.4) X 10"'' 0. The ratio TZf^/^ = 

r(BO^J/^/o(980);/o(980)^.+.-) ^^^^^^^^ j.^ 

< 0.5, based on scal- 



in the range 0.2 < TZ 



fo/'P 

ing from the measurements of Dg decays to /o 
and (f) mesons Using the world-average value 

B{B"g J/V-f/)) = (1.3 ± 0.4 ± 0.2) X 10-3 [4], we obtain 
1.3 X 10-4 < B{B° J/V'/o(980);/o(980) ^ tt+tt-) < 
3.2 X lO-''. A recent study [5| also shows that Tlf^/^ 
can be used to estimate the S-wave pollution in the 
5° — > J/i^cj) analysis of (ig ■ 

We study B^ J/7/'/o(980) in fully reconstructed 



B^ — > J/i/)7r+7r- final states using a 121.4 fb-^ data 
sample collected at the T{hS) resonance with the Belle 
detector at the KEKB collider [Q]. -8° mesons can be 
produced in three T{hS) decays: T(5S') ^ B^B*, B;B° 
and B^B^ where the i?* mesons decay to B^^. The num- 
ber of By'Bl ' pairs in the sample is measured to be 
N^)^(,) = (7.1 ± 1.3) X 10^ using methods described 
in [71, |8[ . Production fractions are measured with fully re- 
constructed — -D77r+ decays as described in [§]. We 
determine the fraction of B*B* pairs among all Bi*^Bi*^ 
events to be fB*B* = (87.0 ± 1.7)% in our full data sam- 
ple. The number of B^ mesons in the dominant B*B* 



B'B' 



production mode is thus iV^o = 2N 
(1.24 ±0.23) X 10^. 

The Belle detector [l^ is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector, 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOP), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect Kl mesons and identify muons (KLM). 



3 



Charged tracks are required to originate within 0.5 cm 
in the radial direction and within 5 cm along the beam 
direction, with respect to the interaction point. Electron 
candidates are identified by combining information from 
the ECL, the CDC {dE/dx), and the ACC. Muon candi- 
dates are identified through the track penetration depth 
and hit patterns in the KLM system. For both electrons 
and muons, the identification efficiency is nearly 100%. 
Identification of charged pions is based on the informa- 
tion from the CDC {dE/dx), the TOE and the ACC. Eor 
a pion from i?° — J/ipfo{980), the momentum- averaged 
identification efficiency is about 96% with a 22% kaon 
misidentification probability. 

Two oppositely charged leptons = e or /i) 

and any bremsstrahlung photons lying within 50 mrad of 
the or tracks are combined to form a J/tp candi- 
date. The invariant mass is required to lie in the ranges 
-0.150 GeV/c^ < Mee(-y) - mj/^ < 0.036 GeV/c^ or 
-0.060 GeV/c2 < Af^^ _ mj/^ < 0.036 GeV/c^, where 
mj/^p denotes the J/^p mass, and Mgg(^-) and are the 
reconstructed invariant masses for e'^e~{"f) and fi'^ ^~ , 
respectively. We combine the J/tp candidate and a tt+tt" 
pair to form a meson. The tt+tt" and J/ip vertex po- 
sitions are required to be consistent. 

Two kinematic variables are computed in the e~^e~ 
collision rest frame: the energy difference AE = Eg — 
£'boamj and the beam-energy constrained mass Mbc — 
V(^bcam)^ ~ (Pb)^i where and p*^ are the energy 
and momentum of the reconstructed -B^ candidates and 
£'beam is the beam energy. To improve the AE and 
Mbc resolutions, mass-constrained kinematic fits are ap- 
plied to J/tp candidates. As B*B* pairs are dominant 
in all Bs Bl events, we focus on the analysis of the 
T(5S') B*B* channel. After an initial loose selection, 
66% of events have multiple candidates. Erom these we 
choose the candidate with an Mbc value closest to the 
nominal B* mass. This requirement has an efficiency of 
90% for the correctly reconstructed signal, according to 
Monte Carlo (MC) simulations. We then select events 
that lie inside a 3tT Mbc signal region for B*B* with the 
criterion 5.4041 GeV/c^ < A/bc < 5.4275 GeV/c^, which 
rejects T(5S') Bs*^ Bs events. We use AE and the 
TT^TT" invariant mass Mt^tt to extract the signal. 

To suppress two-jet-like continuum background arising 
from e+e~ qq {q = u,d,s,c), we require the ratio of 
the second to zeroth Eox- Wolfram moments [llj to be 
less than 0.4. This requirement is optimized by maxi- 
mizing the figure-of- merit Ns / \J Ns + Nb , where A^s is 
the expected number of signal events and Nb is the ex- 
pected number of background events in the (Ai?, Mtttt) 
signal box. Other major background sources are from 
BB {B = B°, B'^, B^) events with one B meson de- 
caying to a final state with a J/tp (denoted J /tpX). We 
use a sample of simulated T(5S') decays, with the most 
recent B meson pair production rates [12l | and all known 



B J /tpX processes, to estimate this background. The 
fit region is chosen to be —0.1 GeV < AE < 0.2 GeV and 
M^^ < 2.0 GeV/c^. Background from B^ J/tPn+n- 
peaks near AE = —0.14 GeV and is outside the fit region. 

A study of the J/ tpX MC simulation is used to cate- 
gorize these background components according to their 
origins and shapes. The expected yields in the fit re- 
gion from each source are: (a) B^ J/tprj', 2.6 events, 
with 77' ^"7 in which the photon is lost; (b) B^ — > 
J/tp{K'^ ,Tr^), 45.3 events, which enter the fit region af- 
ter combining with a random pion; and (c) other J/tpX 
sources, 240.4 events, that do not peak in AE and Mtttt. 
There are negligible correlations between AE and M^ttt 
for (b) and (c), which are parameterized by the product 
of a smooth AE function with a threshold Mtttt func- 
tion for the two-dimensional probability density function 
(PDE). The B+ J/tPK+ and B+ J/tPTT+ back- 
ground shapes are treated separately. For (a), the shape 
and yield are obtained from a dedicated MC simulation 
and the measured branching fraction ^13i], where a MC- 
generated two-dimensional PDE is used, since there are 
correlations between AE and Mt^^ that are difficult to 
parameterize analytically. The non-J/tp background is 
studied with data from a, J/tp mass (M/;(^)) sideband de- 
fined as 2.5 GeV/c^ < Mu^^^ < 3.4 GeV/c^, with the 
regions -0.200 (-0.080) GeV/c^ < M^^) - mj/^ < 



0.048 GeV/c2 for J/tP 



(^+^ ) excluded. The 



shape and yield of the non-J/tp background is obtained 
by fitting and counting the J/tp sideband data. In this 
procedure, the lepton identification requirements are re- 
laxed in the fitting to enhance the statistics. A scale 
factor is used in the counting; this factor is the MC- 
determined ratio of the non-J/tp yield in J/tp selection 
window to the yield in the J/tp sideband region. 

Figure [1] shows the data together with the fitting func- 
tions described below. We find two peaks in the Mt^-k 
spectrum of the events in the AE signal region: one for 
/o(980) and another around 1.4 GeV/c^. We model the 
signal Mt^tt PDF as a coherent sum of a Flatte func- 
tion 'T^ for the /o(980) resonance and a relativistic Breit- 
Wigncr function for a second fx resonance with mass 
moifx) and width Tq^/x): 



7/V 



moifo)^ - M2^ - i{gipi + 9292) 

2 



-ae 



Pj/i,^M^^T{Jx) 



"^o(/ 



X) 



M2 



imo{fx)T{fx) 



, (1) 



where the phase-space factors are pi = 2q/M-^T^^ p2 = 
2(7/f /A/tttt , and the mass-dependent widths arc Fi = 
gipi/m-o and T{fx) = ^Q{fx){q/qQ){mo{fx)/M^^). 
Here q (or go) is the pion momentum in the di-pion rest 
frame where the di-pion mass is M^jr ("^o), while qK is 
the momentum a kaon would have if the resonance de- 
cayed to a kaon pair. The J/tp momentum in the rest 
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frame, pj/-^-, is a phase-space factor outside the modulus, 
and a spin factor inside the modulus for L = 1. The 
Flatte function follows the BES parameterization [isj . 
with the parameters mo(/o) = 965 ± 10 MeV/c^, gi = 
0.165 ± 0.018 GeVVc^ and gaM = 4.21 ± 0.33. The 
AE PDF for signal is parameterized as a sum of two 
Gaussians with width calibrated using a control sample 
of T(5S') ^ BfBf,B° J/'ipK*°{K+TT~) in data. 

Contributions from the self-cross-feed (SCF) events in 
which one or two pion tracks from the signal are mis- 
reconstructed as well as non- resonant 5° — J/i/jtt+tt" 
events are also considered. The SCF component is fixed 
to the MC value of 6.0% of the total signal yield in the fit 
region; the PDF shape is modeled with a non-parametric 
histogram. For the non-resonant background, the Mt^tt 
shape is obtained from a phase-space model and the AE 
shape is the same as that of the J/tpfo signal. 

An unbinned extended maximum likelihood fit to the 
data is performed using the sum of all component PDFs. 
The parameters allowed to vary in the fit are the total 
resonant signal yield, which includes the SCF contribu- 
tion, parameters a, 9, mo{fx), ^oifx), the yield of the 
non-resonant background, and the yield of other J /ijjX 
background. The yield of the J/tlj{K,Tr) background is 
fixed to the MC expectation. 

We obtain 98 ± 15 resonant events corresponding to 
Eq.©, where 63;i^ are from B° 7/^/^/0(980) and 19^^ 
are from B^ — >■ J/ V'/x , and the rest are from the interfer- 
ence. The yield for each signal component is calculated 
using the amplitude squared of that component divided 
by the coherent sum of all amplitudes, with statistical 
errors obtained from error propagation using the covari- 
ant matrix of relevant parameters. The amplitude and 
phase parameters are determined to be a = 0.47 ± 0.10 
and 6 = 1.63 ± 0.98 rad. The fitted mass and width for 
the fx component are mo = 1.405 ± 0.015t°:°J57 GeV/c^ 
and To = 0.054 ± 0.033tg:gj3 GeV, where the first er- 
ror is statistical and the second is systematic. These re- 
sults are consistent with the /o(1370) parameters listed 
in the PDG [3| . Henceforth, we refer to this contribution 
as i?j — J/ipfo{1370); however, the possibility of other 
scalar resonance contributions in this region cannot be 
excluded with the present statistics. The obtained non- 
resonant yield is 4±12, consistent with zero. The yield of 
other J/ip background 262 ±23 is consistent with the MC 
estimate of 240 events. The J/V' helicity distributions, 
shown in Fig. [21 are consistent with a longitudinally po- 
larized J/ip in both /o(980) and /o(1370) signal regions, 
as expected for scalar tttt resonances. 

The signal yields, branching fractions, and signifi- 
cances including systematic uncertainties are listed in 
Table HI The significance is calculated from the log like- 
lihood difference for two parameters in the /o(980) case 
and four parameters in the /o(1370) case, when the cor- 
responding signal amplitude is set to zero. 

Contributions to the systematic error are obtained by 
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FIG. 1: Data fit projections to M„„ for (a) -79.7 MeV < 
AE < -19.7 MeV, and to AE for (b) 0.8 GeV/c^ < A'U^ < 
1.16 GeV/c^ and (c) 1.3 GeV/c^ < M^^ < 1.5 GeV/c^ The 
total PDF is shown with a solid line. The dash-dotted curves 
represent the total background, the dashed curves show other 
J/tp background, and the dotted curves are the non-resonant 
component. 
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FIG. 2: The cosine of the angle 9hci between and the direc- 
tion opposite to that of the in the J/ip rest frame. cos(Shci) 
is projected in the AE signal region and /o(980) and /o(1370) 
signal regions as 0.8 GeV/c^ < M^-^ < 1.16 GeV/c^ (left) 
and 1.3 GeV/c^ < M^r-^ < 1.5 GeV/c^ (right). The expected 
distributions from the fit assuming a longitudinally polarized 
J/ip, which would result from a B° — >■ J/t/j -f scalar decay, are 
superimposed. The curves follow the convention in Fig. [1] 
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TABLE I: Summary of signal yields, significances, and product branching fractions S(-Bs J/ipF-jF tt^tv ), where 
F = /o(980) or /o(1370). 



Mode 




Yield 


Significance 


B{B° J/^F;F^ TV+TV-) 






J/t/./o(980) 


63ti^ 


8.4(7 


(1.16l°:?^(stat.)i°:«,(syst.)«:?|(iV^(.)^(. 


)) X 10-* 




J/V'/o(1370) 


19t^ 


4.2(7 


(o.34i°:}i(stat.)tHi(syst.)t^:S^(A^sJ.)^J* 


)) X 10-* 



TABLE II: Relative systematic errors (in %) for B{B° 
J/tPF;F ^ -K+-K-). 



Source 


i3(F = /o(980)) S(F = /o(1370)) 


AE shape 


1.0 


+0.8,-0.6 


/o(980) shape 


+12.4,-13.6 


+8.9,-5.3 


Background parameters 


+1.9,-1.7 


+1.0,-0.4 


Track reconstruction 


1.3 




Lepton identification 


2.6 




Pion identification 


1.6 




B{J/i> 11) 


0.7 




fB*B* 


2.0 






+22.4, - 


15.5 


Total 


+26.0, -21.1 


+24.5,-16.8 



varying each fixed parameter by its error and are summa- 
rized in Table HH Apart from the Ng,§t normahzation, 
the largest systematic effect arises from the uncertain- 
ties of the Flatte parameters for the /o(980) lineshape, 
where the parameters are varied according to errors in 
[isj . For the signal AE shape, the error on the mean 
value is determined from the beam energy calibrated with 
T{5S) —7- T(lS')7r7r and — ^ Dgn, and the error on the 
width is determined from the control sample. The yields 
of the B'^ — > J/ipT]' and i?+ — > J/'ip{K~^, tt"*") components 
are varied according to the experimental errors on their 
branching fractions. Finally, the non-J/'0 background 
parameters are varied according to the results of the J/ ip 
sideband study. Other background shape uncertainties 
are negligible. 

In summary, we report the first observation of B'^ — > 
J/^/'/o(980) and the first evidence for B'^ J/7/'/o(1370). 
The measured i?" — J/ipfo{980) branching fraction is 
in agreement with the estimate 1.3 x 10"^ ^ S{Bs ^■ 
J/V'/o(980);/o(980) ^ tt+tt-) < 3.2 x lO^^. The signal 
for B° J/V'/o(1370) has a significance of 4.2cr. This 
mode represents a new CP channel that can be used to 
study iJj mixing properties. 

Note added: While preparing the final version of this 
manuscript, we became aware of larXiv:" ! 102.0206, which 
reports similar results for J/V'/o(980). 
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